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THE  INSTABILITY  OF 

LINEAR  HETEROGENEOUS  LANCHESTER  EQUATIONS 


Invariant  properties  of  differential  equations  provide  useful  information 
on  the  nature  of  the  solutions.  The  linear  Lanchester  equations  have  interesting 
invariant  properties  which  this  writer  believes  are  not  readily  known  within  the 
community. 

An  n-dimensional  linear  time  varying  Lanchester  equation  [1]  can  be 
described  by 


X  =  Fit)  x{  t)  ,  xi  t„) 


(1) 


where  F{t)  is  a  nxn  matrix  in  R  and  can  be  represented  by  a  partitioned  matrix 
of  the  form 


F 


\ 


where  the  elements  of  F,  [fj;],  are  all  less  than  1  and  x  is  a  n-dimensional 
"force  strength"  vector  in  R"  and  t,,  is  the  initial  time.  Introducing  the 
fundamental  matrix  [2-8]  with  the  property 


x{t)  =  (j)  it,  tj  X  { tj 

it  can  be  shown  that  the 

t 

f  tr  F{  t)  dt 

det  4)  ( t,  t^)  =  e"°  (2) 


Since  the  trace  of  F,  tr  (F(t)), 


is  zero  for  all  t,  it  follows  that 


det  <])  ( t,  t^)  =  1 


(3) 


for  all  n-dimension  linear  Lanchester  equations. 


Since  in  many  cases  in  practice,  F(t),  is  taken  as  constant  or  piecewise  constant 
over  a  sufficiently  small  intervals  tj  <  t<  t,^.,,  it  follows  from  systems  theory 
that 


4)  (t,  tj  =  T-^  T  (4) 


where  T  is  a  similarity  transformation  and  D  is  a  diagonal  matrix  of  eigenvalues 
for  F(t)  in  t..  <  t  <  t,^,,  assuming  distinct  eigenvalues.  Thus,  it  follows  that 
the 

det  4>  (t,  t^)=det{e^^) 

=  e 


Since  det  ♦(t,  t^)  =  1  we  can  conclude  that 


(5) 


for  all  linear  Lanchester  equations.  It  is  this  property  on  the  sum  of  the 
eigenvalues  which  is  invariant  for  linear  Lanchester  equations.  It  should  be 
noted  that  Eq  (5)  is  also  true  when  the  eigenvalues  are  not  distinct. 

The  analysis  of  Equation  (5)  leads  one  to  conclude  that  either  all 
eigenvalues  reside  on  the  imaginary  axis  resulting  in  harmonic  oscillation  or 
there  exists  at  least  one  eigenvalue  residing  in  the  right  half  plane  so  that  our 
system  is  inherently  unstable. 

When  one  applies  a  finite  difference  approximation  to  the  derivative  in  Eq 
(1)  we  obtain 

X  (t^  +  1)  =  (J  +  F(t^)  )  x(t_i)  (6) 


It  then  follows  that  the  sum  of  the  eigenvalues  (C,),  associated  with  the 
difference  equation,  Eq  (6),  is 


2 


where  n  is  the  dimension  of  x. 

For  Eq  (6)  to  be  stable  all  eigenvalues  must  lie  in  the  unit  circle  so  that 
iCjl  <  1.  otherwise  it  is  unstable.  The  eigenvalues  C,-  associated  with  Eq  (6) 
are  related  to  the  eigenvalues  X-  of  the  continuous  case  by 

Ci  =  l+A, 

SO  that  if  the  continuous  problem  is  nonstable  or  have  all  of  its  eigenvalues  on 
the  imaginary  axis,  then  the  difference  equation  is  inherently  unstable. 

It  is  appropriate  at  this  point  to  present  two  examples  which  demonstrate 
this  invariant  property: 

Case  1:  Linear  Lanchester  equation  of  order  two.  This  is  given  by  the  simple 
equation 


Xi  =  -ax2 

^2  =  - 

(o 

-a"! 

V 

\“ 

P  oj 

which  for  x  =  x^)'^  we  have 

X  =  Fx 

(0  -a\ 


with  F  = 


\ 


0 


/ 


3 


and  a,  ^  real  and  less  than  one. 
The  matrix  F  has  eigenvalues 


=  +  y/a  P 


and 


X_  =  -yja  (5 


so  that  the  2nd  order  system  in  inherently  unstable.  Similarly  the  corresponding 
difference  equation  will  also  be  unstable  since  it  has  an  eigenvalue  of  value 

1  +  /or~P  which  is  outside  the  unit  circle. 

Case  2:  Linear  Lanchester  equation  of  nth  order  (n-50).  This  equation  is  of  the 
form 


X  =  Fx 


where  F  is  50  x  50  and  represented 


by  F  =  "J) 

with  A  and  B  each  being  a  25  x  25  matrix  whose  elements  are  randomly  generated 
by  Matlab  and  less  than  one.  This  means  that  x  is  a  vector  of  dimension  50. 
Using  Matlab  for  the  continuous  case,  it  is  found  that 


50 

A,^  =  1 . 3  X  10‘^^ 

i=l 


with  eigenvalues  in  the  complex  plane  as  shown  in  figure  (1).  It  is  immediately 
seen  that  A,,  have  several  eigenvalues  in  the  right  half  plane  so  that  this  system 
is  inherently  unstable. 


4 


The  finite  difference  equation  for  this  problem  yielded 


V  =1 


50 


consistent  with  our  analysis  and  has  its  eigenvalues  in  the  complex  plane  as 
shown  in  figure  (2).  Please  note  the  scales  in  both  figures  are  different. 
However,  there  are  several  eigenvalues  lying  outside  the  unit  circle  so  that  this 
difference  equation  is  inherently  unstable. 

In  conclusion,  we  have  seen  that  the  structure  of  the  linear  heterogeneous 
Lanchester  equations  leads  to  an  invariant  condition  associated  with  the  sum  of 
the  eigenvalues  with  solutions  which  are  inherently  unstable  or  oscillate 
sinusoidally. 


5 


Figure  2.  Eigenvalues  in  Complex  Plane. 
(Discrete  Case) 


REAL  PART  OF  EIGENVALUE 
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